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We propose a concept of a low-symmetry three-dimensional metamaterial exhibiting a Double- 
Continuum Fano (DCF) optical resonance. Such metamaterial is described as a birefringent medium 
supporting a discrete "dark" electromagnetic state weakly coupled to the continua of two nondegen- 
erate "bright" bands of orthogonal polarizations. It is demonstrated that light propagation through 
such DCF metamaterial can be slowed down over a broad frequency range when the medium param- 
eters (e.g. frequency of the "dark" mode) are adiabatically changed along the optical path. Using 
a specific metamaterial implementation, we demonstrate that the DCF approach to slow light (SL) 
is superior to that of the EIT because it enables spectrally uniform group velocity and transmission 
coefficient. 
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The ability to slow down light to low group velocities 
v g compared with the vacuum light speed c while main- 
taining high coupling efficiency |l| is one of the most 
dramatic manifestations of controlled light manipulation 
in optics. Apart from its fundamental significance, it 
has long-reaching technological applications Q, includ- 
ing enhanced nonlinear effects due to the energy density 
compression by as much as c/v g ; pulse delay and storage 
for optical information processing [3]; optical switching, 
and quantum optics. Most approaches to obtaining SL 
rely on the phenomenon of Electromagnetically Induced 
Transparency (EIT) [4(. EIT and its analogs have been 
demonstrated in several media, including cold [1[, warm 
atomic gases [5[ , and even plasmas [6|, [7| . 

More recently, in response to the emerging applications 
such as bio-sensing, an increasing attention has shifted 
towards obtaining EIT using electromagnetic metama- 
terials [8l-ll0(. Metamaterials enable engineering elec- 
tromagnetic resonances with almost arbitrary frequen- 
cies and spatial symmetries. For example, the EIT phe- 
nomenon has been emulated in metamaterials Ida pos- 
sessing two types of resonances: a "dark" one, which is 
not directly coupled to the incident electromagnetic field, 
and a "bright" one, which is strongly coupled to the in- 
cident field. If the respective frequencies of these reso- 
nances, ujq and gjr, are very close to each other, they can 
become strongly coupled by a slight break of the meta- 
material's symmetry. 

The most serious limitation of the EIT-based SL (EIT- 
SL) stems from the desirability of achieving the SL over 
a broad spectral range, especially for optical buffering 
of ultra-short pulses. Increasing coupling between the 
bright and dark resonances broadens the spectral range, 
yet comes at the expense of increasing both the group ve- 
locity v g and group- velocity dispersion d(v~ 1 )/du. Qual- 
itatively, the bandwidth limitation of the EIT-SL arises 
because the flatness of the EIT transmission band (which 
is necessary for small duj/dk Zl where uo and k z are, respec- 
tively, the frequency and the wavenumber of the propa- 




FIG. 1: (Color online) (a) PBS of DCF metamaterial as 
predicted by the analytical model of coupled oscillators [see 
Eqs. ([Tl[2])]. Different line styles correspond to the reduction 
of the spatial symmetry of a metamaterial: from two mirror- 
symmetry planes (solid) to none (dashed). (b) Broadband SL: 
the medium is comprised of multiple layers of DCF metama- 
terials with spatially-varying resonance frequency Qq in the 
cjo — Acj < Qq < cjo + Acj range. The incoming light under- 
goes polarization transformation and is slowed down. 



gating radiation) originates from the finite spectral width 
2Acj of the EIT band. Therefore, only the light inside 
the ujo — Auj < uj < ujo + Auj bandwidth can be slowed 
down. Using metamaterials with spatially-dependent 
ujo(z) along the propagation direction z cannot increase 
the total propagation bandwidth of the EIT-SL because 
the EIT band is surrounded by the stop bands as illus- 
trated below. Therefore, a natural solution to the band- 
width problem is to design a metamaterial supporting a 
propagating mode which is spectrally broad, yet possesses 
a "flat" segment with du/dk z <C c, as is schematically 
shown in Fig. [TJ 

In this Letter we propose a new approach to producing 
SL which relies on the phenomenon of Double- Continuum 
Fano (DCF) [11] optical resonance: coupling of a single 
discreet state ("dark" mode) to the two sets of contin- 
uum states (two propagating modes of different polar- 
ization states). If the frequency of the "dark" mode is 
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embedded in the frequency continua of the propagat- 
ing modes, then a simple symmetry breaking couples 
all three modes. The result of such coupling is a very 
unusual propagation band [see Fig. [Ha)] which fulfills 
the above requirement for broadband SL applications by 
continuing from one propagating mode to another, with 
the SL region in between. We demonstrate that multi- 
ple layers of such DCF metamaterial with adiabatically 
changing frequency of the "dark" mode [as schematically 
explained in Fig. 0Jb)] give rise to broadband SL with 
spectrally uniform group velocity and transmission. Un- 
like the more commonly known single-continuum Fano 
resonance that can be observed with highly sym- 

metric molecules possessing at least one reflection mirror 
symmetry, DCF resonance requires low spatial symmetry 
molecules that are not reflection-symmetric with respect 
to any plane passing through the direction of light prop- 
agation. 

Before introducing a specific metamaterial realization, 
we consider a very general case of a DCF medium com- 
prised of anisotropic molecules containing two "bright" 
transitions (characterized by the oscillator strengths q x 
and q yi respectively, coupled to the two light polariza- 
tions) and one "dark" transition (characterized by the 
oscillator strength Q and decoupled from the light). The 
equations for the coupled oscillator strengths q x , q y , and 
Q excited by the light with the electric field components 
E x and E y are: 

Qx + QrxQx + Kxy% + KxqQ aa x E x e lUJt : , 



Qy H~ &RyQy H" K xy q x + K y Q( 



aa y E y e iujt , 



Q + VLqQ + K x QQx + l^yQQy = 0, 



(1) 



where £Ir x /r v /q are the resonant frequencies of the tran- 
sitions and K xy / x Q/yQ are the coupling coefficients be- 
tween them. q x and q y are coupled to external fields 
through the coupling constants aa x and aa yi where a 
is a constant with a dimension of cj 2 , and a x / y are the 
depolarization factors such that the dipole moment nor- 
malized to one molecule's volume Vo is p x / y = a x / y q x / y . 
On the other hand, the dark state Q cannot be directly 
excited by the electric field, nor does it directly contribute 
to the polarizability of the medium, thereby playing the 
role of the discrete Fano state. 

Solving for q x and q y in the form of qi = bijEj 
(where i=x,y), we can construct the dielectric permit- 
tivity tensor e = 1 + 47rx, where x = (NoVo)ab, and 
a =diag[<Xc, a y ], where No is the molecular density. The 
effective polarizability of the DCF medium is given by: 
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Xxx = N aa 2 x (W Ry W Q - k\ q 
X V y = N aa 2 y (Wr x Wq - k 2 xQ ) /D 
Xxy = Xyx = -N aa x a y {^ xy W Q - k x qk v q) /D 



(2) 



where D = W Rx WR y WQ-Kl Q W Ry -nl Q WR X -Kl y WQ- 
^KxQKyQKxy is the determinant of the matrix formed by 
the l.h.s of Eq. JT), and W { = (ft 2 - a; 2 ). Photonic band 
structure (PBS) is analytically calculated from the eigen- 
value equation for k z (uj): det(k z 5ij—6ijUJ /c 2 )=0, which 
directly follows from the Maxwell's equations. 

Unlike the EIT case where the dark mode's frequency 
needs to be matched to that of the x-dipole resonance [8|, 
[Tot , the DCF occurs under the following conditions: (i) 
Qq < ^R X ^R y: and (ii) at least 2 of the 3 coupling 
constants (frxQiKyQ, an d K X y) are non- vanishing. With- 
out significant loss of generality, and with an eye on the 
specific metamaterial implementation (see schematic in 
Fig. [2]), we have neglected the direct coupling between q x 
and Q. Finite values of the ft's is the consequence of the 
reduced spatial symmetry of the anisotropic molecule. 
For example, n xy = would be either accidental, or due 
to the high spatial symmetry of the molecule (e.g., y — z 
mirror symmetry). Figure [TJa) shows the PBS in the 
DCF medium (dashed line) calculated from the analytical 
model using the following parameters: [Qr x , ^R y , ^q] = 
[27,23.85,19.1] THz, [^,^,^ = [50,0,50] THz 2 , 
a x =l, a y =1.65, a=529 THz 2 , and N V =10- 3 . 

The unusual PBS in the DCF medium can be under- 
stood through its evolution from that in the medium 
comprised of the symmetric molecules (k x q = 0, n xy = 0) 
[solid lines in Fig. QJa)]. First, the k x q ^ coupling 
hybridizes the x-dipole and the quadrupole resonances 
and creates an avoided crossing between the x-polarized 
propagation band and the discrete "dark" mode. The 
second avoided crossing between the ^/-polarized propa- 
gation band and the discrete "dark" mode occurs owing 

tO K X y ^ 0. 

Thus, the DCF medium comprised of the low- 
symmetry molecules supports a propagation band that 
smoothly connects the two orthogonally polarized con- 
tinua by passing through the SL region near the "dark" 
mode's frequency Qq. Such PBS has the following ad- 
vantages over that in the EIT medium: (i) the dark 
mode frequency has more tunability because it is only 
required to be embedded inside the continuous bands, 
and (ii) there are no bandgaps on either side of the 
SL region. These features of the SL band are neces- 
sary for realizing the broadband SL using the approach 
of adiabatically- varying material parameters [specifically, 
Qq(z)] schematically shown in Fig. QJb). 

Both DCF and EIT metamaterials supporting the 
propagation of slow light can be implemented using a 
three-dimensional periodic arrangement of the unit cells 
comprised of three metallic antennas. Such unit cell 
shown in Fig. [2] consists of 2 vertical and 1 horizontal 
metallic antennas embedded in a dielectric medium with 
n = 1.5. It is reminiscent of the ones used for single-layer 
EIT metamaterials [ilflQj]. except that both s x and s y are 
allowed to be non- vanishing, thereby dispensing with all 
reflection symmetries of the unit cell. The two vertical 
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caused by the band- folding owing to the periodic nature 
of the metamaterial in the z-direction. 
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FIG. 2: (Color online) PBS for SL metamaterials based on 
the DCF resonance (a) and on the EIT (b). The insets 
show the geometry and dimensions of the unit cell and the 
three supported resonances: (i) horizontal dipole, (ii) vertical 
dipole, and (iii) the quadrupole. Solid lines: PBS computed 
for a symmetric unit cell (s x = and s y = 0). (a) Prop- 
agation bands for DCF-based metamaterial (dashed lines): 
s x = 700 nm, s y = 2 /xm, L± = 2 /mi. The avoided crossings 
marked (1) and (2) are caused by s y / and s x / 0, respec- 
tively Flat portion of the spectrally-extended middle band: 
SL. (b) Propagation bands for the EIT-based metamaterials 
(dashed lines) with partial symmetry breaking (s x — nm, 
s y = 500 nm, L± = 3.8 //m): emergence of SL for the 
spectrally-narrow middle band. For (a) and (b): t =400 nm, 
h =400 nm, L 2 =4 /mi, w =800 nm, g 2.2 /im. The 
metamaterial's periodicities are 6 /mix 7 /mix 7 /mi, and the 
electromagnetic waves are assumed to propagate in the z di- 
rection. 



antennas support a bright dipole and a dark quadrupole 
resonances, while the single horizontal antenna supports 
another bright resonance. Within the general model de- 
scribed by Eqs. (pQ), the strengths of these resonances 
are measured by q y , Q, and q x , respectively. Coupling 
between all 3 resonances is induced by breaking the re- 
flection symmetries of the structure. For example, s y ^ 
results in k x q ^ 0, while s x ^ results m hv X y 

The PBS, calculated using the finite-elements software 
COMSOL, is shown in Fig. [2] for different geometric pa- 
rameters of the unit cell. For example, by approximately 
matching Qr x and Qq (this is done by choosing the ap- 
propriate length Li of the horizontal antenna) and se- 
lecting s x = 0, s y 7^ 0, an EIT-SL band is obtained 
as indicated by the middle dashed line in Fig. EJb) (y- 
polarized bands are not shown). Note that the EIT-SL 
band is surrounded by the two stop-bands. On the other 
hand, when (i) s x ^ 0, and (ii) the dark mode of the 
fully-symmetric structure (s x = s y = 0) intersects the 
two propagation continua (solid lines), a DCF-SL band 
emerges as shown in Fig. [2fa) (dashed lines). The dif- 
ference between the DCF bands in Figs.QJa) and[2fa) is 
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FIG. 3: (Color online) Field enhancements of four differ- 
ent frequency components propagating in an adiabatically 
varying DCF-based metamaterial. Inset: single- layer unit 
cell. Spectral position of the SL is adiabatically varied: 
L2[/im] = 3.75 + z/700. Field intensity is calculated at the 
red spot shown in the inset. Other parameters: s y = L2/4, 
w — t — 0.8 /xm, g — 2.7 /mi, and gap — 1.7 /mi. 

While it is apparent from the PBS that light can be 
significantly delayed by either EIT or DCF-based meta- 
material, the spectral width of such SL band is quite 
limited. To create a broadband SL, we propose a multi- 
layered structure with each layer having its DCF's res- 
onance frequency ~ Qq adiabatically varied along the 
propagation direction z. In the metamaterial structure 
shown in Fig. [2] such tuning can be achieved, for exam- 
ple, by changing the length of the 2 parallel antennas. 
Because the SL mode in the DCF-based metamaterial 
contains no band gaps, the adiabatic process smoothly 
converts one fast mode (e.g., x-polarized) into the hybrid 
slow mode, and then into another fast mode (^/-polarized) 
without significant reflection. Such conversion is not pos- 
sible for EIT-SL because the SL band is surrounded by 
the stop bands which lead to reflection of the SL as was 
recently demonstrated 3, 15]. 

Specifically, we consider a broadband x-polarized laser 
pulse with the central frequency c^o and bandwidth Auj 
incident upon an TV-layer metamaterial with adiabati- 
cally varying Qq^ < Qq^ < fig^ (where j is the meta- 
material's layer number) as shown in Fig. QJb). It is 
further assumed that, while Auu is much larger than the 
spectral width of the SL portion of any individual layer, 
the following relations are satisfied: < ujq — Auj and 

coo + Auj < . As the pulse passes through the struc- 
ture, each frequency component uj is slowed down inside 

(i) 

its corresponding layer j satisfying Qq w uj while propa- 
gating with no significant delay through the other layers. 
The light pulse is also gradually converted from x- into 
^/-polarized. Because all the frequency components of the 
pulse undergo the same adiabatic transition, the entire 
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pulse is slowed down uniformly with no group velocity 
dispersion. 

COMSOL simulations of a TV = 41-layer DCF-SL 
metamaterial, with the single-layer unit cell shown in the 
inset to Fig. [3j were performed. Note that this unit cell, 
while topologically equivalent to the one shown in Fig. [2j 
is optimized to decrease Ohmic losses. Antennas are as- 
sumed to be made of silver with the dielectric permit- 
tivity described by the Drude model [16| . The structure 
is triply-periodic with the period L = 7 jim and embed- 
ded in a dielectric with n = 1.5. The SL's frequency is 
adiabatically varied from layer to layer by increasing L2 
from 3.75 jam to 4.15 jam. Other unit cell parameters are 
given in the caption. Field intensity enhancements eval- 
uated at the corner of the left vertical antenna are shown 
in Fig. [3] for 3 different frequencies. As Fig. [3] shows, 
different frequencies are slowed down inside different lay- 
ers as indicating the corresponding field enhancement. 
Thus, light is slowed down over the entire spectral band 

ng>< w <n™ 
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FIG. 4: (Color online) Transmission (solid green lines), ab- 
sorption (dashed red lines), and group velocity (solid blue 
lines) of (a) adiabatic DCF-based and (b) EIT-based meta- 
materials. x-polarized incident light is assumed. Parame- 
ters of the DCF-based metamaterial: the same as in Fig. [3] 
EIT-based metamaterial is made of 20 identical layers shown 
in Fig. [2] with the parameters L2 = 4 /mi, L\ — 3.8 /mi, 
w — t — 0.8 /xm, g — 2.7 /xm, h — 0.9 /mi, s y — 0.6 /xm, and 
periodicity L x = L y — L z = 7 /mi.. 

To underscore the advantages of the DCF-based ap- 
proach to broadband SL, we compare a multilayer DCF- 
SL structure (unit cell is shown in Fig. [3]) with its EIT- 
based counterpart (unit cell shown in Fig. [2j geometric 
parameters are given in the caption to Fig. 3]). Both 
structures are based on silver antennas embedded in an 
n = 1.5 dielectric. Figure |4] (a) and (b) shows the COM- 
SOL simulations results for DCF-SL and EIT-SL struc- 
tures, respectively. The group velocity is calculated as 
Ldu/d(j), where L is the total length of the structure and 
(j) is the phase difference between incident and transmit- 
ted waves. In the case of EIT, the bandwidth can be 
potentially increased by increasing the coupling between 
dipole mode and quadrupole mode (k> x q). However, it 



comes at the expense of the increased group velocity dis- 
persion and transmission non-uniformity. 

Specifically, both v g and transmittance |T| 2 increase in 
the middle of the EIT band while remaining small near 
the EIT band edge. In contrast, the adiabatic DCF layers 
provide uniform group velocity as well as uniform cross- 
polarized transmittance inside the SL band. The band- 
width can be increased further simply by adding more 
adiabatically- varying layers. Note that the group veloc- 
ity in the DCF-SL structure plotted in Fig. [4] is averaged 
over N = 41 layers. Inside the specific resonant layer, the 
group velocity of the corresponding frequency component 
is < 0.01c, which is consistent with the intensity enhance- 
ment of three orders of magnitude. The spectrally-flat 
transmission and absorption of the DCF-based metama- 
terial is appealing to a variety of linear and nonlinear 
applications. Moreover, because every frequency compo- 
nent is dramatically slowed down in a well-defined layer, 
applications to spectrally-selective active light manipula- 
tion can be envisioned. 

In conclusion, we have proposed a new mechanism 
of slowing down light over a spectrally broad band by 
means of low-symmetry metamaterials exhibiting double- 
continuum Fano (DCF) resonance. This approach is con- 
ceptually different from the more common EIT-SL. It is 
shown that DCF-based broadband slow light with uni- 
form group velocity and transmittance can be achieved 
by adiabatically changing the DCF resonance frequen- 
cies along the light propagation direction. Spectrally-flat 
light absorbers and filters, as well as various nonlinear 
devices requiring extreme light concentration are enabled 
by DCF-based slow-light metamaterials. 

This work was supported by the grants from the Air 
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search. 
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